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Abstract 

We compute the annihilation cross-section of two neutralinos at rest into two photons, which 
is of importance for the indirect detection of neutralino dark matter in the galactic halo through 
a quasi-monochromatic gamma-ray line. We include all diagrams to one-loop level in the minimal 
supersymmetric extension of the standard model. We use the helicity formalism, the background- 
field gauge, and an efficient loop-integral reduction method. We confirm the result recently 
obtained by Bergstrom and Ullio in a different gauge, which disagrees with other published 
calculations. 
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1 Introduction 


The composition of dark matter is one of the major issues in astrophysics. A popular candidate for 
non-baryonic dark matter is a neutral stable Majorana fermion, the lightest neutralino, appearing in 
a large class of supersymmetric extensions of the Standard Model. In a wide range of supersymmetric 
parameter space, relic neutralinos from the Big Bang are in principle abundant enough to account 
for the dark matter in our galactic halo. (See ref. Q for a thorough analysis and ref. for a 
comprehensive review.) 

Several methods have been proposed to detect galactic neutralinos: elastic scattering in a low- 
background terrestrial detector, energetic neutrino fluxes from neutralino annihilations in the core 
of the Sun or of the Earth, gamma-rays and cosmic ray antiprotons or positrons from neutralino 
annihilation in the halo. (For a review see ref. ||^].) 

In particular, the annihilation into two photons gives a distinct signature against backgrounds 
from known astrophysical sources - a narrow line in the gamma-ray spectrum at the energy equal 
to the neutralino mass |3|. Observation of such a line would provide evidence for the existence of 
supersymmetric dark matter, while non-observation can be used to put constraints on models of the 
galactic dark matter. A number of new experiments with reasonable energy resolutions are being 
planned Q which would be sensitive to photon production in the galactic halo. 

There exist several previous calculations of the neutralino annihilation into two photons, all but 
one incomplete. The first calculation was carried out by Bergstrom and Snellman for neutralinos 
which are pure photinos. Rudaz corrected their results and considered also pure higgsinos. Giu- 
dice and Griest |Q removed the restriction on the neutralino state mixture, and confirmed Rudaz’s 
work. In all these early calculations the sfermion and the VE-boson masses were assumed to be 
much larger than any other mass. This assumption was relaxed by Bergstrom Q, who computed 
the fermion-sfermion loop contributions for arbitrary squark masses, but only for a pure photino. 
Bergstrom and Kaplan evaluated the contribution from virtual IE’s in a leading-logarithm ap¬ 
proximation. Jungman and Kamionkowski improved on the leading-logarithm approximation, 
but did not include all contributing diagrams (our below was missing). Only very recently, 

Bergstrom and Ullio |l^ (hereafter BU) presented a full one-loop calculation. Their results, however, 
differ from some previous partial calculations when corresponding expressions are compared. 


We provide an independent complete one-loop calculation of the neutralino annihilation into 
two photons. We confirm the BU result using a different calculational procedure. To perform the 
calculation we take advantage of some recent improvements in the calculational techniques for one- 
loop amplitudes |n]. In particular, we use a helicity basis \12, for the photons, background field 
Feynman-’t Hooft gauge 15] and an efficient integral reduction method 


As in previous calculations, we focus only on the first term in the relativistic expansion of av, 
where a is the annihilation cross-section and v the relative neutralino-neutralino velocity. This 
kinematic limit is appropriate for neutralino annihilation in the galactic halo, where the neutralino 
velocities are of the order of the galactic rotation speed, u/c ~ 10“^. This kinematic configuration 
considerably simplifies the form of the amplitudes. Our evaluation methods are however general and 
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do not depend on the special kinematics. 


2 Evaluation of the Amplitudes 

In this section we outline the computation of the amplitudes for the annihilation of two neutral 
massive Majorana fermions at rest into two photons via vector boson, scalar and fermion loops. The 
Feynman diagrams for the process under consideration are given in fig. 0. Since photons do not 
change the identity of the particle they couple to, the masses in the fermion and scalar lines within 
the loops are uniform. 



Figure 1: Feynman diagrams for neutralino annihilation into two photons. 


2.1 Feynman Rules and Diagrams 


We use background field Feynman-’t Hooft gauge [14, because of its technical advantages. In 
this gauge the gauge fixing part of the Lagrangian is 


^gf = - 


1 

2 





1 

2 






2 


( 1 ) 


where W and B are respectively the SU{2) and 17(1) external background fields, W and B are the 
corresponding quantum fields appearing in the loops, and (pok is the VEV of the k-th. Higgs field of 
the model: (pk = 4>ok + (p'k- Fadeev-Popov ghosts are rather simple and at one-loop have the 
same coupling as scalars; only the overall sign differs because of Fermi statistics. The Feynman rules 
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Figure 2: Feynman rules in the background field gauge; overall factors of i and coupling constants have 
been removed. 


for performing this calculation may be found in ref. |^] and are summarized in fig. In this figure 
7 and Z are background fields. 

An important technical advantage of this gauge is that only one of the three terms in the WW^ 
vertex contains loop momentum. This decreases the number of the integrals with high powers of 
loop momentum in the numerator that have to be evaluated. Typically, such integrals are the most 
difficult to compute. Furthermore, the gauge fixing terms have been chosen so that they cancel 
the WH'j coupling from the Lagrangian, which reduces the number of diagrams to be evaluated. 
(This is similar to non-linear gauges 0, Si .) This cancellation works for any number of 
Higgs fields assuming that they are doublets of the Standard Model SU{2). For the massive four- 
point amplitudes computed in this paper, the background field method captures the main technical 
advantages of using the gauge invariant cutting techniques reviewed in ref. |11|. 


2.2 Helicity and Kinematics 

We use a helicity representation for the final-state photons, since this generally leads to compact 
expressions for amplitudes. In the spinor helicity formalism [^, |l^ the photon polarization vectors 
are expressed in terms of massless Weyl spinors |A:^) = ^(l± 75 )|/c). 


'ik-,q) = 


(g 

\/2 (qk) 


:{k]q) = 


y/2 [k q] 


( 2 ) 


where k is the photon momentum, q is an arbitrary null ‘reference momentum’ which drops out of 
final gauge-invariant amplitudes and 


{ij) = (A;. |A;+), 


[ij] = {kt\k- 


(3) 


which satisfy the normalization condition {ij) [ji] = 2ki ■ kj. The plus and minus labels on the 
polarization vectors refer to the outgoing photon helicities. We denote the massive Majorana spinors 
of the initial-state neutralinos by | 1) and | 2). (For the massive spinors we choose the normalizations 


of ref.[181; our massless spinors follow the conventions of ref. |13||.) 
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Since the two annihilating particles are identical Majorana fermions, the initial-state wave- 
function must be antisymmetric. Therefore the initial spin state must be the singlet with angular 
momentum zero; since it is conserved, the final-state photons must also have vanishing angular mo¬ 
mentum and be of the same helicity. This means that there are only two helicity cases that are 
required, namely where both photons are of positive or both are of negative helicity. These two 
helicity configurations are, however, simply related as can easily be seen by making the reference 
momentum choice 

q3 = h, qi = ks, (4) 


where and are the reference momenta of the two photons. With this choice of reference momenta 


t/l — 


(3 4) 




— 


(34) 




[34] 3 ’ ^3 [ 34 ] ^4 • 

Since the amplitudes are symmetric under the interchange of photons 3 4, we have 

^X’ ^7 ^^' y ) — r .,2 "^(^X’ ^X’ 3 ^) , 


[3 4]^ 


(5) 

( 6 ) 


where legs 1 and 2 are the initial state neutralinos and the ± labels on legs 3 and 4 refer to the 
photon helicities. Thus, the amplitudes for either helicity choice are identical up to an overall phase 
factor and we need only compute the annihilation into two positive helicity photons to obtain the 
entire result. 


A number of important simplifications occur because of the special kinematics. For annihilation 
at rest, we take the kinematic conhguration to be of the form 


/ci = (mj^, 0,0,0), A:2 = (m^, 0,0,0), 

fes = (-m^, 0,0, m^), /c4 = 0, 0, -m^), 

where is the neutralino mass, fei and /c 2 are the neutralino four momenta and and k^ the 
photon momenta. The Mandelstam variables are then given by 


s = (fci-I-A: 2 )^ = 4m|, t = {ki + k^)'^ =, u = {ki + k^)'^ =. ( 8 ) 


The four-point kinematics allows us to reduce any spinor structure that appears in the diagrams 
(before we symmetrize in the photons) to a linear combination of (1|2), ( 1 | 75 | 2 ), (l||b 4 | 2 ), and 
(1| 75^4 |2). The permutation 3 4 then corresponds to k^ ^ k^, t ^ u. (Note that the reference 

momentum choice (Q) respects the exchange symmetry.) But in our kinematics t = u, so we only 
need to symmetrize the spinors. Using (1| + ^4 |2) = 0, (1| 75 ( 1(3 -|- ^ 4 ) |2) = 2my^ (1| 75 |2), we 
see that the answer after symmetrizing in the photons contains only (1|2) and (II 75 I 2 ). At the 
kinematic point (^) we have ( 1 | 2 ) = 0 , so all amplitudes are proportional to ( 1 | 75 | 2 ). 

Another simplification due to this kinematics is that we need only consider pseudo-scalar Higgs 
particles and Z bosons in the s-channel, the other neutralino currents vanish. This means that 
in fact we do not need to calculate diagrams and in the diagrams we need to consider 

only pseudo-scalar intermediate Higgs. Moreover, diagrams A'^^ do not contribute because CP 
conservation (which we impose) forbids the and couplings. 
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2.3 Integral Reductions 


In evaluating the diagrams in fig. one encounters loop integrals of the form, 


In[P{n] = -^(4vr) 


2-e 


( 27 r) 4 - 2 e LILI..LI 


where 


T 2 /?2 2 

L, = £• — m, . 


( 9 ) 


( 10 ) 


Here the ii and rrii are the momentum and the mass flowing through the rth loop propagator, and 
P{i^) is a polynomial in the loop momenta. For the box (four-point) integrals, 

£i = i , £2 = £ — ki, £3 = £ — ki — k2 , £4 = £-|-A:4 . (11) 

(Our sign conventions for the integrals In are adjusted to agree with Oldenborgh’s FF program, and 
differ from those of refs. 0, 0, 0] for the three-point integrals.) Although, the amplitudes are 
ultraviolet finite we include dimensional regularization because a few of the integrals encountered in 
the calculation are divergent. (We comment that we use the ‘four-dimensional helicity scheme’ |21] 
to maintain compatibility with the helicity formalism; however, since the amplitudes are ultraviolet 
finite the specific scheme choice is unimportant.) 

We now outline the procedure used to reduce the tensor integrals (i.e. those with powers of 
loop momentum in the numerator of eq. (^)) to linear combinations of scalar integral functions 
used to express the amplitudes. This procedure has been used previously in refs. |I^, |2^, |ll|. The 
basic technique is to extract as many inverse propagators as possible from the spinor inner products 
appearing in the numerators of the integrands. In this way we can simplify the tensor box integrals 
with two or three powers of loop momenta, which are by far the most complicated of the integrals 
occurring in the calculation. 

As an example, consider the box diagram in fig. |T](a) with a scalar cj) and a charged fermion ■i/’ in 
the loop, given by 




( 1 | ^2 + + i^\2) ■ {£4^ + £ 3 ) ef- (£i + £4) 


J (27r)4 (27r)-2^ {ij - - ml){£l - iP - ml){£l - ^P - rin?^){£l - iP - m^) ' 

As usual when using four-dimensional helicities p^ , 22, 0, we have broken up the loop integral into 
a four-dimensional part, and a (—2e)-dimensional part, See, for example, appendices A.2 and 
C of ref. 122] for further details. 

Using the polarization vectors d), with the choice of reference momentum (Q) the dot products 
of the polarization vectors and loop momenta may be re-expressed as 

(3+K4I4+) (4+1413+) 


63 • {£4 + £3) 64 • {£i + £4) — —2 ■ 


(3 4)^ 


(13) 


We may then extract inverse propagators from the spinor products using, 
(3+|/4|4+)(4+|/4|3+) 

= ( 3 +I/ 4 K 4 / 4 I 3 +) 

= 24 • 4 (3+ 1/4 |3+) - ^ 4 ( 3+1 ^4 |3+) 

= [(^4 -ml- ^P) - (£f -ml- ^P)]2 4 • ^3 - [{£l -ml- fp) + ml + fP]s34 


(14) 
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where the combinations appearing in the parentheses are inverses of propagators appearing in the 
loop integral (12). By canceling propagators we can reduce the tensor box integral (12) to a linear 
combination of box integrals with at most one power of loop momentum and triangle integrals with 
at most two powers of loop momentum. In some cases we can reduce the tensor triangle integrals by 
again extracting inverse propagators; however, this strategy fails for some terms because one cannot 
form an inverse propagator appropriate for the triangle integral. For example, in eq. (0) , the term 


'I-ml- 1X^)2 ■ ks = {il -ml- [{^l -ml- i?) - {il -ml- /i^)] 


(15) 


contains two factors of (^ 4 —/i^) which cannot be completely canceled against a single propagator. 
For such terms one can use, for example, the integration formulas of refs. [^, |T^, ^ to directly 
evaluate the tensor integrals in terms of scalar integrals; the expressions are relatively simple since 
these are triangle integrals. 

The integrals with a single power of ^ in the numerator vanish, while the ones containing a 
are proportional to e times a higher-dimensional integral [p 2 |, pi]] , 


(27r)4 (27r)-2^ 


y? = -dvre 


(27r)6 (27r)-2^ 




(16) 


The higher dimensional integral may be re-expressed in terms of H = 4 — 2e integrals using the 
integral relations appearing in refs. |W, j^. Observe that because of the explicit e, we need only 
evaluate the parts of the higher dimensional integrals that are singular in e. In some cases, such 
as for box integrals containing a single power of the integrals are finite and may therefore be 
dropped immediately since they are suppressed by an overall power of e. It turns out that for the 
amplitudes presented in this paper, all such higher dimension integrals cancel. 

The above procedure allows us to reduce all four-point diagrams in fig. Q to a linear combination 
of scalar box, triangle and bubble integrals. 


2.4 The Cross-Section 


We calculate the cross-section in the framework of the minimal supersymmetric extension of the 
Standard Model. (For a comprehensive discussion, see ref. |^].) We use the following notation for 
the MSSM interaction terms that contribute to the diagrams in fig. |I|: 

Ant = + + ^-c. 

+ + h.c. (17) 

+ + h.c. 

Here y, are Majorana fermions, are Dirac fermions, and cf) can be a neutral {(j)* = (j)) or a 

charged scalar. Explicit expressions for the coupling constants gijk can be obtained in refs. |^ , |^] . 
We write the amplitude as 


A{1^,2^,2,+ AX) 


4iam|(l|75 |2) - 

O 

IT (34)^ 


(18) 
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where a is the fine structure constant. With this notation, A matches the corresponding expression 
in refs. |9|, 

The cross-section may be directly obtained from the amplitudes listed in eq. (^) below using 
the formula 

a^m% ~ „ 

= 

We have averaged over the four possible initial states, summed over the two possible photon helicity 
configurations and inserted the symmetry factor of ^ for identical final state photons. In evaluating 
the squared matrix elements we have used 


1(34)1 = 2m^, |(1|75|2)| =2v/2m;^, 

valid for the kinematics (D- 

We find the following non-vanishing contributions to the amplitude; 

-^= E E + E+ E E -4: 

'>P=Xi,xt <t>=H±,G± <p=f (t>=HlG° ■>P=f,xt,xt 


( 20 ) 


Zip 


( 21 ) 


'tp=xt,xt 


^=f iJ=f,xt^xt 


where the sums are over all the Standard Model fermions / (quarks and leptons), the corresponding 
sfermions /, the charginos xf, and the unphysical Higgs bosons G^, with masses rriQO = mz 
and mQ± = mw, which appear in’t Hooft-Feynman gauge. The individual amplitudes are 

= “2 Swxp} 2 4^^ {mw) + 2{ml + mw - m|) 4“' -h 2m^ 4^' + 4""' + 4^^ i^w, m^) 

+ 8 Dwx'>p rn^m^{lf^ -h 4""') > 

2m^ 4^' + "^0 4""' + 4^^ ^i>) + 4^' > 


0 V' 

id 


Ad Z r r 




4m? 

X 


m-^ m?^ — 4m? ^ 


4^^(m^) 


2 A A ^0 4m? 

•^10 = C 0 e ^ 52 xx fi '^00 4^2 _ ^2 h ("^ 0 )' 


( 22 ) 

where is the electric charge of particle xjj in units of the positron charge, is a color factor (3 
for quarks; 1 for leptons and charginos). 


Swjk = + \gfjk\'^ , Dwjk = \gijk\'^ - \g, 


12 


- IW |2 


12 

ijk\ 


S<pjk = \gfjk? + \g[jk\‘^, D^jk = - h 


r.P |2 


S |2 




(23) 


The labels for amplitudes follow the labels appearing in fig. |I|. For example, is the gauge 

invariant combination of diagrams appearing in fig. |^(a), (b) and (c) which have a W boson in the 
loop. Similarly, is the same combination of diagrams, but with the W in the loop replaced 

by a scalar. 

For the triangle diagrams, A'^^r and do not contribute because of the special kinematics 
(see discussion in sect. ^)^) . A^w^ “4^^ and the part of A^-^ with the vector coupling between Z 
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and the fermion loop vanish by Furry’s theorem. The term requires special attention because 
of the anomaly, which arises when the Z coupling with the fermion loop is axial. We have dropped 
the anomalous contributions in eqs. (|^ ) since they add to zero when summing over all fermions in 
the loop. The remaining contribution of this diagram is non-zero due to the large mass splitting 
between the two quarks of the third generation and between the charginos. 

Individual Feynman diagrams are, of course, not gauge-invariant. However, the combinations 
appearing in eq. (|2^ are gauge invariant. We have verified that when the photon polarizations are 
taken to be longitudinal these combinations vanish as required by gauge invariance. This provides 
a strong check on our calculation, including relative factors between diagrams. 

Substituting the at rest kinematics (^) into the general integral identity given in eq. (18) of 
ref. we can reduce the box integrals appearing in to triangle integrals: 


r[&] 


2 I 9 2 ’ 

2(2 2 ’ 
'^X + 


(24) 


The integral functions in are obtained from eqs. (^) by replacing mw with 

We use the standard methods of ref. to express the triangle integrals in terms of logarithms 
and dilogarithms |28], 


/W(m) = < 




1 -|- X\1 


8 m| - 

X 


2m| ( 


m < 


arctan 


/rn^ — 1 / 


4 ^'(mi, m 2 ) = 


2 m| 


Li 2 (^ 

Li 2 (| 


mf — m| — m^ — 


2 m? 


) + Li2(^ 


mf-I-m| - m| - ■y/Ai's ^ / 
2mf ^ 


m > m^; 

mf — m| — mf -|- \ 

2mf ^ 

mf m| - mf -h \/A^ \ 

2mf ^ 


where 


(25) 


(26) 


X = \J'^ — m?/rnj^, 

Ai = (mf -k m? — mf)^ -k 4m?mf, 

A 2 = (mf — m| — mf)^ — 4m|mf. 

One convenient way to obtain the values of the integrals is from Oldenborgh’s program FF [^]. In 
the notation used by FF the triangle integrals are 


l 3 ^^(m) = /sfm^, m^, m^; 4m?, 0,0), 

4 ^'(mi,m 2 ) = -m?,0 ,m?). 


(27) 
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We have compared the results in eq. (^) to those of BU and we find complete agreement. We 
also find disagreement between our results and those of ref. |p. After using eq. ( p4|) to express the 
amplitudes in terms of triangle integrals, the comparison is simple. For example, for the integrals 
of class 1 and 2 of BU the translation between the real parts of our integrals and those of BU 


is. Re 13 ^^( 777 .^) /i(a,6)/4m|, — 13 ( 0 , 6)/2m| and ^ —l2{a,b)/2m 

(Notice that 13 ( 0 / 6 , 1 / 6 ) = 12 ( 0 , 6 ).) Similar formulas hold for the class 3 and 4 integrals. The 
notation for the coupling constants is related by g^Pi + Q^Pr = + < 7^75 (for vector bosons; 

g^ + g^lb for scalar bosons) with Pl = {1 — 75)/2 and Pr = (1 + 75 )/ 2 . (Note that BU use unitary 
gauge for the Z in the diagrams in fig. ||(e), while we use Feynman-‘t Hooft gauge.) 


r[2] 




4- 


3 Discussion 

3.1 Large mass behavior of amplitudes 


An interesting feature of our result is that in the limit when the neutralino is a pure higgsino with 
a mass much greater than mw the cross-section does not fall off when is increased, but remains 
constant. This behavior, which was numerically observed by BU, is caused by the diagrams in 
%• ^(b). In the pure higgsino limit, x is almost degenerate with the lightest chargino and the 
denominator of in eq. ( p^ ) becomes small, so this contribution dominates. This allows us to 
obtain an analytic expression in the limit ~ S> mw- Using eqs. (p4D, (25), and the 

expansion 

(28) 


Li 2 (l -I- ie) + Li 2 (l - fe) ~ —— evr, 

O 


the asymptotic behavior of is 


74 — 


2 rri'^mw 

This yields the asymptotic behavior of the cross section. 


(29) 


4 

a TT 


av 


10"^®cmVs. 


4m^ sin^ 9w 


(30) 


where we have used = < 7^/4 = vro:/ sin^ 0 ^ and set the other couplings to zero 

which is appropriate for pure higgsinos. 


3.2 Annihilation to gluons 

We may also use the above results to obtain the cross section for neutralino annihilation into two 
gluons. The diagrams for this process are a subset of the ones appearing in the photon calculation. 
(The diagrams containing three-gluon vertices vanish.) This allows one to obtain the result for gluon 
annihilation directly from the results for photon annihilation. For the gluon case, the amplitude is 
given by 

■^ = T ^57""’ + U ■iii. + Y. -^7. (31) 

,l>=q <t>=H°,G° 1 /’='? 

i/i=g 
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where we set = 1 in the individual amplitudes and sum over quarks and squarks. The color 

sum amounts to simply changing c? to in (p!^), so our result agrees with BU. This confirms a 
small but significant error in previous calculations |]30| , , which generally led to an over-estimation 

of antiproton fluxes^ 


4 Summary and Conclusions 


In this paper we presented a complete one-loop calculation for neutralino annihilation into two 
photons or two gluons within the Minimal Supersymmetric Standard Model. In performing the 
calculation we assumed that initial-state neutralinos are non-relativistic, which is the relevant case 
for galactic dark matter. The techniques we used are, however, not limited to this special kinematics. 


Our results confirm the very recent calculation of Bergstrom and Ullio [IT] which disagrees with 
previously published ones, especially in the case where the neutralino is a heavy pure higgsino. 


We thank Graciela Gelmini for encouragement and a number of stimulating discussions. This 
work was supported in part by the DOE under contract DE-FG03-91ER40662 and by the Alfred 
P. Sloan Foundation under grant BR-3222. P.G. thanks Piero Ullio and Lars Bergstom for useful 
discussions and David Cline, Graciela Gelmini and Roberto Peccei for the extended visit at UCLA 
which made this work possible. 
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